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ABSTRACT 

We report a new solution to explain the observed deficit of the solar neutrino 
flux by Homestake, Kamiokande II and III, GALLEX and SAGE experiments. We 
use the matter mixing and the helicity oscillation in the twisting magnetic fields 
in the sun. Our model predicts the short (seasonal) and long (11 years) time 
variations of the solar neutrino flux. Three kinds of data observed by Homestake, 
Kamiokande, GALLEX and SAGE detectors are reproduced well if the mixing 
angle and the squared mass difference are in the small area around sin 2 26 ~ 0.01 
and Am 2 ~ 1.3 x 10~ 8 eV 2 . 



1. Introduction 



Three kinds of data of the solar neutrinos flux are now available from Homes- 
takefl], Kamiokande II and III[2], GALLEX[3] and SAGE[4] experiments. These 
data have shown the deficit of neutrino flux with respect to the prediction of the 
standard solar models [5], [6], [7]. The different deficit rates with respect to detec- 
tors show the neutrino energy dependence of the observed flux. If we take all 
of these data seriously, there needs new physics. The MSW (Mikheyev-Smirnov- 
Wolfenstein) matter mixing is the most natural scenario [8] of new physics. 

Another characteristic feature is that the Homestake experiment seem to see 
the anti-correlation between the number of sunspots and the counting rates. This 
anti-correlation was not observed by the Kamiokande II and III experiments. Var- 
ious authors [9] have examined this and concluded that the anti-correlation was not 
significant. Recently, we proposed a new interpretation of the time profiles of the 
Homestake and Kamiokande II data[10]. We take the complex time dependence of 
the Homestake data as a real phenomena and consider that the data show both 
the short (seasonal) and long (11 years) time variations. We interpreted the short 
time variation is originated from the effect of the twist of toroidal magnetic fields 
and the long time variation (11 years) is due to the change of their magnitudes. 
The Kamiokande II data are the averaged ones for about one year so that no short 
time variation shows up. 

We constructed a simple model of the twisting toroidal magnetic fields following 
the simulation by Yoshimura[ll] and made the numerical analysis. We successfully 
reproduced the Homestake data as well as the Kamiokande II data[10]. Since we 
are intended to see purely the effect of the twisting toroidal magnetic fields, we did 
not take into account of the matter mixing effect. Thus our predictions did not 
have the neutrino energy dependence. As a result, the GALLEX and SAGE data 
are not be explained. 

In this paper, we consider the model which have both mechanisms; the matter 
mixing and the helicity oscillation in twisting toroidal magnetic fields. Our concern 
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is to find a new kind of solution which explain both the deficit of the neutrino flux 
(the time averaged profile and the neutrino energy dependence of the flux) and 
the short and long time variations of the flux. We made the numerical analysis 
and examined the allowed region in (Am 2 , sin 2 26) plain. We found that if Am 2 ~ 
10~ 8 eV and sin 2 26 ~ 0.01, all data are reproduced well. 

In Sec. 2, we discuss the evolution equation and its general features. A model of 
the twisting toroidal magnetic fields are briefly explained in Sec. 3 and the general 
features of oscillation of this scheme is discussed in Sec. 4. Numerical analysis and 
the allowed region are given in Sec. 5. Sec. 6 is devoted to the discussions. 



2. The evolution equation 

Our model consists of two Majorana neutrinos, u e and which have the tran- 
sition magnetic moment fi. In the following, we use v e and for left-handed 
neutrinos and v e and for their anti-neutrinos, respectively. The evolution equa- 
tion for neutrinos which fly along z axis under magnetic fields is given by 



.dip 
dz 



( V e 

S 2 5 





s 2 5 

-V^ + 2c 2 5 fiBre* 
nBre-** -V e 
s 2 5 





S 2 5 

V fI + 2c 2 5j 



4, 



where i\) T = (v e , v^, v e , u^), s 2 = sin 26 and c 2 = COS26 1 with 6 being a mixing angle, 
5 = Am 2 /AE with E being the energy of neutrino and Am 2 being the squared mass 
difference. The magnetic fields relevant to the helicity oscillation is the transverse 
component which are parametrized by B^e 1 ^ = B x + iB y . If the twist of toroidal 
magnetic fields exists, has the z dependence. The matter potentials for v e and 
are denoted by V e = Gp(2n e — n n )/y/2 and V^ = Gpn n /\/2 where Gp is the 
Fermi constant, n e and n n are number densities of electrons and neutrons in the 



3 



sun, respectively. For them, we take [12] 



V e =O.195i?0 1 exp[lO.82 v / l - z/Rq] 
=0.018 j Rq 1 exp[10.82 v / l - z/R Q }. 

This is the extension of the OVV (Okun, Voloshin,Vysotsky) model[13]. The 
model without the matter mixing (0 = case) was discussed by various au- 
thors^], [15], [10] by considering various types of twisting magnetic fields. The 
application of this model to the solar neutrino problem with a realistic model of 
twisting magnetic fields was made in Ref.[10]. 

The importance of the phase can be understood by changing phases [14], 
ip T — > ip T = (e*^/ 2 z/ e , e^/ 2 z/ M , e~^/ 2 z/ e , e~*^/ 2 z/ M ). Then, apart from the overall 
phase, ip obeys the following equation 



dz 



( Ve s 2 5 -fiB T \ 

s 2 S -V^ + 2c 2 5 -ilB t 

(j,B T -V e -<j)' s 2 5 

\-fj.B T s 2 5 Vf, + 2c 2 5 -0' J 



(3) 



where 0' = d(f>/dz. The important role of the phase lies in the fact that the 
variation of it works as a potential. Due to this, various kinds of resonance oscil- 
lations will occur as discussed by Akhmedov, Petcov and Smirnov[16], depending 
on sizes of parameters, c 2 5, s 2 5, /iBt and 0'. There are essentially two resonance 
oscillations of v e . In the following, we assume Am 2 > and cos 29 > because we 
are interested in the transition of v e . 

(i) The MSW matter (flavor) oscillation 

If fiBx is small in comparison with s 2 5, the MSW matter oscillation occurs. 
The resonance condition for v e <-> oscillation is 

V e + Vf, = (Am 2 /2E) cos 29. (4) 
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When the adiabatic condition at the resonance point 



1 



d(V e + 



dz 



Am 2 sin 2 29 
2E cos 29 



(5) 



is satisfied, v e converts fully to v^. If this condition is not respected, the conversion 
is not full and the transition rate is estimated by the numerical computation or 
the Landau-Zener formula. The v e <-> oscillation does not occur because the 
condition V e + V^ = — (Am 2 /2E) cos 29 is not met. 

(ii) The OVV helicity oscillation in the twisting magnetic fields 

If S2# = (Am 2 /AE) sin 20 is smaller than /i-By, the OVV oscillation in the 
twisting toroidal magnetic fields occurs. In this case, the oscillation between two 
different flavor neutrinos occurs because Majorana neutrinos can have only transi- 
tion moments. The resonance point is determined by 



Ve-Vfj, = (Am 2 /2E)cos29 - <p' (v e <-> P^), 



(6) 



\/ e -V^ = -(Am 2 /2£)cos2#-0 / (p c <-> 

We consider the survival probability P of v e after passing the distance Az. The 
rate is expressed by [17] 



P(v e -> u e ; Az) = 1- B f sin s 



^(V e - V„ - 2c 2 S + <P>) 2 + WBt) 2 



Az, (7) 



where Bf is the blocking factor defined by 
Bf 



4(/x5 T ) 2 



(8) 



(V e - V^-2c 2 5 + <P') 2 + A(iiB T ) 2 - 

The above formula is valid only when the variation of V e — V^ + <p' is small in 
comparison with /i£>y. In this formula, the matter potential acts as a blocking 
factor. At the resonance point, the blocking factor disappears (Bf = 1) and the 
conversion occurs. The adiabatic condition 



d(V e -V^ + 00 



dz 



< A(nB T f 



(9) 



at the resonance point. 
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Since the energies of neutrinos spread from about 0.33MeV (pp neutrino) to 
about lOMeV ( 8 B neutrino), 5 = Am 2 /4i? varies about factor 30 so that neutri- 
nos may receive different types of effects. We look for the situation where the pp 
neutrinos receive mainly the MSW matter oscillation effect, while the 8 B neutri- 
nos receive the OVV helicity oscillation effect in the twisting toroidal magnetic 
fields. In this situation, the time variations occur for 8 B neutrinos, but not for pp 
neutrinos. This possibility is realized when 5 and fiB^ are the same size. 

3. A model of twisting toroidal magnetic fields 

Here we give a brief summary of the model of twisting toroidal magnetic fields 
which is proposed in our previous papers[10]. As in Fig.2 in Ref.10, We assume 
two tori in the convective area of the sun; one in the northern hemisphere and the 
other in the southern hemisphere parallel to the equator. We assume that toroidal 
magnetic fields locate in the torus and toroidal magnetic fields twist along it. This 
twist is parametrized by a parameter X, the distance along the torus to wind once. 
As in Fig. 5 in Ref.10, we parametrize the configuration of the torus in the sun's 
cross section. The torus is parametrized by the latitude of its center A, its radius 
a, the distance between the center of torus and the center of the sun b. 

The latitude of the neutrino path is parametrized by A which is between —7.25° 
(the southern hemisphere) and 7.25° (the northern hemisphere). Since A is small, 
neutrinos pass through almost around the edge of the toroidal magnetic fields so 
that we assume that the strength of magnetic fields is constant along the neutrino 
path. We specify the position of the entrance zo of a neutrino to magnetic fields 
and the departure z\ from it. They are given by 

20,i = b cos(A - A) ± sja 2 -b 2 sin 2 (A - A). (10) 

The twist of troidal magnetic fields is generated by the differential rotation of 
the sun and the global convection of the plasma fluid in the convective area. By 
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solving the dynamo equation, Yoshimura[ll] showed that the development of the 
twist of the troidal magnetic fields is the origin of the cyclic oscillations of polarity 
reversals in every llyr. The direction of the twist depends on whether magnetic 
fields lie in the northern or the southern hemisphere. The global structure of the 
twist is determined by the Coliori's force which acts on the global convection. By 
comparing the simulation by Yoshimura with our simple model, we found that the 
variation of the phase <$' is given by [10] 

// / \ \ 2ir/X 

{Z) = - Sign(A) l + ( 27 r/X)2[&cos(A-A)-,]2 - < U ) 

The characteristic feature is that < if neutrinos pass through the northern 
hemisphere and (f)' > if they do the southern hemisphere. This sign difference 
will give the important affect on the deficit rate of the neutrino flux[10]. Various 
parameters defining the toroidal magnetic fields are fixed by comparing our model 
with the simulation by Yoshimura, A = 15°, a = O.1694i? and b = O.8813i? . 



4. General features of oscillation 

We see some qualitative features of the neutrino survival rate. In order to make 
some concrete arguments, we take sin 29 = 0.1, (/i-Br) ~ 4 x 10~ 10 /iekG~ 8/Rq 
and X = tiRq (about two turns) so that <$' ~ ±2/i?©. Also, we restrict A in two 
cases A = 7° (northern hemisphere, around September) and A = —7° (southern 
hemisphere, around April). With these values, we obtain zq = 0. 758/2© and z\ = 
0.99LR©. 

The survival rate is given as a function of y which is defined as 

(E/lMeV) 3.5 x 10 9 E 
V = (Am 2 /leV 2 ) ~ Rq Am 2 /(leV 2 )' ^ 

The evolution equation is solved from the center to the surface of the sun with the 
initial condition ip = (1,0,0,0) at the center of the sun where v e is created. In 
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Figure 1, we show our result about the y dependence of the survival probability 
P{y e — > u e ) observed on earth. The solid (dotted) line represents the survival 
probability of v e when it passes through the northern (southern) hemisphere. The 
seasonal difference arises for y > 2 x 10 7 where the effect of the twisting magnetic 
fields becomes dominant. In Figure 2, we show to which spices u e transformed. 
Fig.2b(c) corresponds to the case where neutrinos pass through northern (southern) 
hemisphere. v e transforms mainly to up to y ~ 3 x 10 7 . When y is greater than 
this value, u e transforms to v^. The results in these figures may be understood 
qualitatively as follows: 

(i) y < 2 x 10 5 

The resonance condition of the flavor oscillation in Eq.(4) is not satisfied and 
the resonance oscillation does not occur. 

(ii) 2 x 10 5 < y < 10 6 

In this region, the resonance point is inside the sun and the adiabatic condition 
is satisfied. Thus, v e transforms fully to v^. The resonance condition of helicity 
oscillation may be satisfied, but the resonance point is no in magnetic fields which 
are located in the region 0.758-R© < z < 0.991-R©. Thus, the helicity conversion 
does not occur. 

(iii) 10 6 < y < 3 x 10 7 

In this region, the non-adiabatic flavor transition occurs. As in Figs, lb and lc, 
a part of v e transforms into v^. 

(iv) 3 x 10 7 < y < 5 x 10 8 

Since /jBt is assumed to be around S/Rq, it is larger than (Am 2 /4i?) sin2^. 
Thus the helicity transition v e — > occurs. The resonance condition is expressed 
by Ye = 0.177 exp(10.82A/l — zJRq) and the adiabatic condition in Eq.(9) is 
expressed by 



where Ye — 1.8 x 10 9 /y — <J)'Rq. In order for the adiabatic condition is satisfied, 
| Y# | <d 1 is required. This means that the resonance point should be very close to 
the surface and thus the adiabatic condition is not satisfied. In this non-adiabatic 
case, the estimation of the transition probability is rather complicated. We see 
from Fig. la that in this y region the seasonal difference appears due to the effect 
of the twist 4>' . The role of (f>' is to shift the resonance point. When neutrinos 
pass through the northern (southern) hemisphere, (f>' is negative (positive) and the 
resonance point z r is smaller (larger) than the non-twisting case. For example, 
for y ~ 10 8 , the resonance point determined from Eq.(7) is z T ~ 0.81i? Q for 
the northern hemisphere and z x ~ 0.837?© for the southern hemisphere. Since the 
transition is non-adiabatic, the transition develops gradually until neutrino reaches 
to the surface. It turns out that for y ~ 10 8 and (f>' < (northern hemisphere), 
the survival rate becomes minimum before neutrino reaches to the surface so it 
oscillates back at the surface. On the other hand, for 0' > (southern hemisphere) 
the minimum value is achieved at the surface. This is the reason why the seasonal 
difference appears. 

(v) 10 9 > y > 5 x 10 8 

In this region, the pure helicity transition in the matter discussed in Ref.10 is 
realized. Mostly, v e transforms to v^. The general tendency is that the suppression 
of v e flux is strengthen when neutrinos pass through the northern hemisphere 
((f)' < 0) and is weaken when neutrinos pass through the southern hemisphere 
((p' > 0). The transition rate depends sensitively on the strength of magnetic fields 
as shown in Ref.10. 

Since the average energy of pp neutrinos is 0.33MeV and that of 8 B neutrinos is 
lOMeV, there is about 30 times difference. The experimental data by Homestake, 
Kamiokande II and III, SAGE and GALLEX suggest that pp neutrinos are less 
suppressed than 7 Be and 8 B neutrinos. From Fig. la, there seem two possible 
solutions: One is the region y ~ (10 6 ~ 10 7 ) where the MSW matter oscillation 
works. The pp neutrinos are less suppressed than 7 Be and 8 B neutrinos and no 
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time dependence of fluxes appears. In this solution, v e transforms to v^. The other 
is the region y ~ (10 8 ~ 10 9 ) where the pp neutrinos are less suppressed and have 
small time dependence. On the other hand, the 7 Be and 8 B neutrinos are more 
suppressed and two kinds of time dependencies, long (llyr) and shot (seasonal) 
oscillation appear. This is the region where we are seeking. 



5. A new solution 

In order to solve the evolution equation, we have to give the size of the twist cf>' 
and the time variation of the strength of magnetic fields. Although the twist should 
exist, the size of <p' is an unknown parameter. Here we take X = ttRq. We guess 
the time variation of magnetic fields from the variation of the sun spot number, 
Bt oc V sunspotnumber. We also assumed the maximum strength (/x.B:r) max = 
10.5 x 10 _10 /iBkG. With these choices of values, we can estimate the neutrino 
flux at each year and each season. For the latitude of the neutrino path, we 
consider A = 7° (around September) and A = — 7° (around April). According to 
the simulation by Yoshimura, magnetic fields on neutrino paths from August to 
October are similar to those of September and thus the calculation with A = 7° 
will valid for the period from August to October. Similarly, the calculation with 
A = — 7° will be valid for the period from March to May. 

We now seek a solution in the region y ~ 10 8 following the discussion given in 
the previous section. We considered three experiments, Homestake, Kamiokande 
II and III and GALLEX. For Kamiokande II and III, we used the formula 

P(KII) = P( 8 B) + 0.11P( 8 B), (14) 

where P and P represent the probabilities of v e and 9^, respectively. The symbol 
8 B in the parenthesis show the survival rate of 8 B neutrinos which is evaluated with 
use of the average energy E u =10MeV. The contribution from the anti-neutrino 
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arises because Kamiokande II and III detectors are sensitive to it also. GALLEX 
and Homestake detectors can catch other than 8 B neutrinos. We used the formula 



P(GALLEX) 



70.8 




132 



(15) 



P (Homestake) 



0.2 
7^9 




(16) 



For pp, 13 N and 5 neutrinos, the survival rates are estimated by using their 
average energies 0.33, 1.1 and 1.44MeV. 

For Kamiokande II and III and GALLEX data, we demand that the time aver- 
ages of our calculated rates agree with their averages within 1 standard deviation. 
We relaxed this constraint for the Homestake data because the data points changes 
with time so that the time average may not give a good measure for agreement. 
We put the importance on the time variation more than the average value for this 
case. We require that the time average of our estimates agree with the average 
experimental value within four standard deviation. 

Allowed regions are shown in Figure 3. The area covered by dash-doted lines, 
solid lines and dotted lined are the allowed regions by GALLEX data, Homestake 
data and Kamiokande II and III data, respectively. We found an unique shaded 
area which satisfies three experimental restrictions. The mixing angles and the 
squared mass in the allowed domain are sin 2 26 ~ 1 x 10~ 2 and Am 2 ~ 1.3 x 



We show the comparison between our predictions and the data in Figure 4. 
Fig. 4a show the comparison with Homestake data, Fig. 4b with Kamiokande II and 
III data and Fig. 4c with GALLEX data. Our predictions for the Homestake data 
show the short (seasonal) and long (llyr) time variations and reproduce the data 
pretty well. We also predict these time variations for the Kamiokande II and III 



10~ 8 eV 2 . 
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data, but the comparison with respect to the short time variation will be premature 
at present because the Kamiokande II and III are given as the average for about one 
year. There appears no seasonal time variation if one year average data are taken. 
As for the GALLEX data, the main component is pp neutrinos. Thus we predict 
almost no time variation for GALLEX data as seen in Fig. 4c. The agreements of 
our predictions and three kinds of data seem to be pretty good. 



6. Discussions 



We presented a new solution to explain the solar neutrino problems, the missing 
flux problem and the short and long time variations. We used a simple model of 
twisting magnetic fields in the sun and showed that this type of magnetic fields 
can give a new kind of solution for the solar neutrino problem. Our model predicts 
the short and long time variations which should appear in the Homestake and 
Kamiokande II and III data, and weak time dependence for GALLEX and SAGE 
data. Our predictions seems to reproduce these data pretty well. We also estimated 
to which type of neutrinos u e converts. From Figs. lb and lc, the pp neutrino 
converts mainly to v^, while the 7 Be and 8 B neutrinos do to v^. Thus, the situation 
in our model is more complicated than the MSW matter oscillation case. In order 
to clarify these situations, we have to wait the future experiments. 
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FIGURE CAPTIONS 



1) Energy dependence of the survival probability of v e . We estimated it for 
sin 26 = 0.1, (fJ.B T ) ~ 4x 10 -1 VflkG an( j x = vri?©. The solid (dashed) line 
corresponds to the case where neutrinos pass through the northern (southern) 
hemisphere. The MSW matter oscillation mechanism dominates when y = 
(E /lMeV) / (Am 2 /leV 2 ) is smaller than 3x 10 7 so that no seasonal difference 
arises. When y > 3 x 10 7 , seasonal differences appear due to the twist of 
toroidal magnetic fields 

2) Energy dependence of neutrino components. We show where v e transforms 
when neutrinos pass through the northern hemisphere (a) and the southern 
hemisphere (b). For y < 3 x 10 7 , v e transforms mainly to v^, while for 
y > 3 x 10 7 , v e does to v^. 

3) Allowed domains by three experiments. Homestake data restrict the area to 
the domains surrounded by the solid lines; a large area in the upper-right 
corner and many small islands. Kamiokande II and III data restrict to the 
areas surrounded by the dotted lines; relatively large areas with Am 2 ~ 
l(T 8 eV 2 and ~ (2 ~ 3)l(T 9 eV 2 and several islands. GALLEX data restrict 
to the large area from Am 2 ~ 3 x 1CP7 to 4 x 10 8 eV 2 which is surrounded 
by the dash-dotted lines. The intersection of these areas is the small area 
with the shade. 

4) The comparison of our predictions with three data. With the choice of some 
values of Am 2 and sin 29 in the shaded area, we estimated the neutrino yields 
for three experiments. The comparison of our prediction with Homestake 
data is shown in (a), with Kamiokande II and III data in (b) and with 
GALLEX data in (c). 
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